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Abstract. There is increasingly strong observational evidence that slow magnetoacoustic modes 
arise in the solar atmosphere. Solar magneto-seismology is a novel tool to derive otherwise 
, directly un-measurable properties of the solar atmosphere when magnetohydrodynamic (MHD) 
wave theory is compared to wave observations. Here, MHD wave theory is further developed 
illustrating how information about the magnetic and density structure along coronal loops can 
be determined by measuring the frequencies of the slow MHD oscillations. The application to 
observations of slow magnetoacoustic waves in coronal loops is discused. 
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Damped slow MHD oscillations have been observed in the solar atmosphere using 
high-resolution EUV imager onboard space-borne telescopes (see review by Wang 2011). 
Such oscillations are important because of their potential for the diagnostics of magnetic 
' structures by implementation of the method of magneto-seismology, through matching 
J> , the MHD wave theory and wave observations in the solar atmosphere to obtain several 
physical parameters (e.g., magnetic field strength and density scale height). 

The theory of MHD wave propagation in solar magnetic structures initially began 
modelling the magnetic structures as homogenous cylindrical magnetic flux tubes en- 
closed within a magnetic environment (Roberts et al. 1984). Later on, more advanced 
equilibrium models to study slow MHD oscillations have also been proposed with, e.g., 
dissipative effects and gravity (Mendoza-Briceno et al. 2004, Sigalotti et al. 2007), and 
non- isothermal profiles (Erdelyi et al. 2008), while the effect of density and magnetic 
stratification had been revisited on transversal coronal loop oscillations by Dymova & 
Ruderman (2006) and Verth & Erdelyi (2008), respectively. 

Here, the governing equation of the longitudinal mode is solved numerically for density 
stratified loops with uniform magnetic field, as well as for expanding magnetic flux tubes 
with uniform density. The effect of these stratifications on the frequency ratio of the first 
overtone to the fundamental mode is studied. 



2. Governing equation 

The ideal MHD equations are linearized by considering small magnetic and velocity 
perturbations about a plasma in static equilibrium [b = (b r ,0,b z ) and v = (v r ,0,v z ), 
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for r and z the radial and longitudinal coordinates, repectively] . In the derivation, a 
uniform kinetic plasma pressure is assumed, and the thin flux tube approximation is 
considered. The second-order ordinary differential equation governing the longitudinal 
velocity amplitude is (see Luna-Cardozo et al. 2012 for a detailed derivation) 
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where c\ = (B 2 /pp ), c 2 = {jpo/po), c\ = c 2 s +c\ and c\ = (c s 
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[ 2 ) 1 are the square 

of the Alfven, sound, fast phase and tube speeds, respectively. In this equation, uj is the 
angular frequency of the oscillations. 

Equation (|2.1[) is numerically solved using the shooting method based on the Runge- 
Kutta technique, for density stratified loops with uniform magnetic field, as well as for 
expanding loops with uniform density. Solar waveguides are modelled as axisymmetric 
cylindrical magnetic tubes with tube ends frozen in a dense photospheric plasma at 
z = ±L. On average, plasma density and magnetic field strength are expected to decrease 
with height above the photosphere (Lin et al. 2004). 



3. Effect of density stratification 

The solar coronal loop is modelled by a straight axisymmetric magnetic flux tube with 
tube length of 2L and radius of ro- The uniform magnetic field is directed along the tube 
axis, i.e., B = B z z. In semi-circular coronal loops where the plasma is close to hydrostatic 
equilibrium, a reasonable assumption for the density profile is 



p Q (z) = pi exp 
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where H is the density scale height and pf the density at the footpoint. To study a 
standing wave the boundary condition v z (±L) — is applied. We solve equation (I2.ip 
using the density profile (|3.1[) . The frequency ratio of the first overtone to the fundamental 
mode is shown in Figure [IJa) as a function of L/H by the solid line, and it is clearly 
lower than the cannonical value of two. A similar result was obtained for the transversal 
mode by Dymova & Ruderman (2006) and Verth (2007). 

For vertical chromospheric flux tubes the density profile is given by 



Po{z) = Pf exp 



(z + L) 



H 



(3.2) 



Longitudinal oscillations in chromospheric flux tubes are studied solving the eigenvalue 
problem in half of the magnetic bottle, i.e., designating v z (—L) — v z (0) = as the 
boundary conditions. The ratio of frequencies against L/H for the density profile (|3.2|) is 
shown by the dashed line in FigureUJa). Now, the frequency ratio is slightly greater than 
two, indicating that this parameter depends on the functional form chosen of the equilib- 
rium density. This suggests that caution must be used when interpreting the frequency 
ratio of chromospheric standing modes. 



4. Effect of a non-uniform magnetic field 

An expanding flux tube with rotational symmetry about the z-axis in cylindrical coor- 
dinates (r, 9, z) is used to model a magnetic field equilibrium decreasing in strength with 
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Figure 1. (a) Frequency ratio of the first overtone and fundamental mode against L/H for 
density stratified coronal (solid line) and chromospheric (dashed line) loops, (b) Frequency ratio 
against the expansion parameter V for different values of fit in vertical chromospheric flux tubes. 
In (b) dotted, dot-dashed, dot-dot-dot-dashed and long-dashed lines correspond to fit — 1, 2, 5 
and 10, respectively. 

height above the photosphere. The magnetic field component B z at the tube boundary 
can be described explicitly as function of z (see Verth & Erdelyi 2008) 



B z {z) « B z>i \ 1 



(1 - T 2 ) [cosh^/L) - cosh(l)] 



(4.1) 



T 2 1 - cosh(l) 

where T — r a /r{ is the expansion factor, and r a (rf) is the apex (footpoint) radius. The 
loop expansion has been estimated for various loops, giving mean values of T » 1.16 and 
1.30 for EUV and soft X-ray loops (Watko & Klimchuk 2000, Klimchuk 2000). 

We can compute the numerical solution of equation (|2.ip for slow longitudinal oscil- 
lations in coronal and chromospheric loops setting the same boundary conditions as in 
the previous section, and using equation (|4.ip for B z (z). FigureQJb) shows the frequency 
ratio as function of the expansion parameter T for different values of the footpoint beta 
plasma /3f for chromospheric flux tubes. It is found that when the magnetic expansion 
factor increases the frequency ratio decreases, and this effect is more significant for chor- 
mospheric flux tubes with higher /3f . 

Figure [2] shows the frequency ratio as function of the expansion factor for coronal 
loops with uniform density in (a) and for typical density stratification (i.e., L/H = 2) in 
(b). It can be seen how these two effects, density stratification and magnetic expansion, 
contribute to decrease the value of uj^/i^x- Additionally, the effect of the expansion is 
stronger in the corona than in the chromosphere. 



5. Summary and conclusions 

Studying the solutions of the velocity governing equation of the slow standing mode, 
it is found that density stratification and magnetic expansion cause the same qualitative 
effect on the frequency ratio in coronal loops, giving values of < 2. For chromo- 

spheric flux tubes density stratification and magnetic expansion cause opposite effects on 
the frequency ratio; however, caution must be taken when studying chromospheric flux 
tubes since the ratio a^/wi depends on the functional form chosen for the density (see 
Luna-Cardozo et al. 2012 for an analytical and numerical detailed study about important 
issues). 
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Figure 2. Frequency ratio of coronal loop oscillations against the expansion parameter F for 
different values of /3f . Solid, dashed, dotted, dot-dashed and dot-dot-dot-dashed lines correspond 
to fit = 0.1, 0.5, 1, 2 and 5, respectively. Coronal loops with uniform density (L/H = 0) are 
presented in (a) and with typical density stratification (L/H = 2) in (b). 

These results are consistent with the values of period ratio of P1/P2 = 1.54 and 1.84 
reported by Srivastava & Dwivedi (2010) while observing slow acoustic oscillations using 
Hinode, in contrast to the theoretical value of P1/P2 = 2 for a uniform cylindrical flux 
tube model. 

Our results arc important for magneto-seismology, where the density scale height of 
the solar atmosphere can be calculated by using the observed value of the frequency 
ratio uj 2 /uji of longitudinal loop oscillations, to complement both emision measure and 
magnetic field extrapolation studies. This could provide us with a more complete un- 
derstanding of the plasma fine structure in the solar atmosphere. These results can be 
applied in any stage of the solar cycle, including the solar minimum. 
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